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Abstract 

The Tolman-Bondi solution of the Einstein equations is used in order to 
model the time evolution of the void observed in Bootes. The present den- 
sity contrast of the central region (~ —0.75) and its radius (~ 30/1" 1 Mpc) are 
fixed, while the density parameter of the Universe, the amplitude of the density 
contrast inside the void wall, the width of this wall and the distance from the 
void centre to the Local Group are appropriately varied. The microwave back- 
ground anisotropy produced by Bootes-like voids is estimated for a significant 
set of locations. All the voids are placed far from the last scattering surface. 
It is shown that the anisotropy generated by these voids strongly depends on 
the density parameter, the wall structure and the void location. The Doppler 
dipole and quadrupole are subtracted and the residual anisotropy is calculated. 



In the case of some isolated Bootes-like voids placed at redshifts between 1 and 
10 in an open universe with density parameter Slo = 0.2, the residual anisotropy 
appears to be a few times 10~ 6 on scales of a few degrees. This anisotropy is 
about one order of magnitude greater than previous estimates corresponding to 
other cases. The anisotropy produced by a distribution of voids is qualitatively 
studied in the light of this result. Comparisons with previous estimates are 
discussed. 
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1 Introduction 

Two methods have been used in order to estimate the anisotropies produced by the 
nonlinear voids of the galaxy distribution in open universes. One of these methods is 
based on the so-called Swiss-Cheese model (Rees & Sciama 1968). The original version 
of this model applies to the case of overdensities surrounded by underdensities, but 
suitable modifications lead to a model for underdensities surrounded by overdensities 
(Thompson & Vishniac 1987). Estimates of the anisotropies produced by Swiss- 
Cheese voids were obtained by Thompson & Vishniac (1987), and Martinez- Gonzalez 
& Sanz (1990). The main elements of these spherical voids are: an absolute vacuum 
in the void core, a uniform overdense shell compensating the vacuum, and a general 
Friedmann-Robertson- Walker background outside this shell; hence, the density profile 
of the resulting structure is very particular, while the background is general and the 
compensation is exact. The second method is based on the Tolman-Bondi Solution 
(TBS) of the Einstein equations (Tolman 1934; Bondi 1947). A general assymptotic 
Friedmann-Robertson- Walker background and a general spherically symmetric energy 
density profile are compatible with this solution; such a general profile can be used 
in order to model both a partial vacuum in the void core and a sharped wall. Two 
different codes based on the TBS were built up by Panek (1992) and Arnau et al. 
(1993). Panek's code was used (Panek 1992) to estimate the anisotropy produced 
by voids with small compensating walls, while the code due to Arnau et al. was 
used (Arnau et al. 1993) in the case of voids without walls; here, the void walls are 
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modeled in some detail taking into account some observational data. 

In flat universes, another powerful method is being used in order to estimate 
the anisotropy produced by nonlinear cosmological structures. This method does 
not require any symmetry, it is based on the potential approximation developed by 
Martinez-Gonzalez, Sanz & Silk (1990), which applies beyond the linear regime. In the 
flat case, the method was applied by Anninos et al. (1991), Martinez-Gonzalez, Sanz 
& Silk (1992, 1994), Tuluie & Laguna (1995) and Quilis, Ibahez & Saez (1995); these 
authors used various complementary techniques and conditions (N-body simulations, 
high resolution shock capturing methods and particular spectra and statistics). 

Recently, Arnau, Fullana & Saez (1994) proved that some Great Attractor-like 
structures -evolving in open universes (Qo < 0.4) and placed at redshifts between 
2 and 30- produce anisotropies of the order of 10 -5 on angular scales of a few de- 
grees. These structures had a density contrast of the order of 1 when they influenced 
the microwave photons; hence, the resulting nonlinear gravitational anisotropies are 
produced in the mildly nonlinear regime. The following question arises: Are there 
void-like objects -suitable structures, locations and backgrounds- producing signifi- 
cant anisotropies as in the case of the Great Attractor-like objects studied by Arnau, 
Fullana & Saez (1994)?. In the case of nonlinear structures (density contrasts 5 > 0.1) 
evolving in open universes, the TBS and the Swiss-Cheese model can be used in or- 
der to answer this question. The use of the TBS seems to be preferable because this 
solution involves appropriate density profiles. 
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In order to model void-like objects, some observational data must be taken into 
account. In the eighties, there has been a great deal of observations on the spatial dis- 
tribution of galaxies; in these observations, some voids with walls have been detected 
(Kirshner et al. 1981, Davis et al. 1982, Vettolani et al. 1985, de Lapparent, Geller 
& Huchra 1986, Rood 1988, Dey, Strauss & Huchra 1990); among them, the Bootes 
Void (BV) seems to be the greatest one. This void was first described by Kirshner 
et al. (1981). From the data given by these authors and those due to de Lapparent, 
Geller & Huchra (1986) it follows that the void in Bootes is a big quasispherical re- 
gion with a defect of galaxies; its centre is located at ~ 150 Yr 1 Mpc from the Local 
Group and its radius is ~ 30 hr l Mpc. This region only contains the 25% of the 
galaxies expected in the same volume of the background; hence, the energy density 
contrast of galaxies inside the void is ~ —0.75. Surrounding this region, there is an 
irregular shell having an excess of galaxies, in other words, there is an inhomogeneous 
sharped wall. Although the observations are not very accurate (Kirshner et al. 1981, 
de Lapparent, Geller & Huchra 1986), current data suggest that the amplitude of the 
density contrast of galaxies inside the wall is ~ 4 and the mean width of this wall is 
~ hh- 1 Mpc. 

The above observational data must be complemented with appropriate assump- 
tions about the dark matter distribution. Since the nonlinear anisotropy produced 
by voids located far from the last scattering surface is a gravitational effect, this 
anisotropy is produced by the total energy density contrast — . This contrast can be 
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obtained from the observational value of the density contrast produced by galaxies 
(^f-)gai an d the value of the the so-called linear bias parameter b. This parameter is 
defined by the relation {^ £ -) ga i = K^r)- ^ n this paper, it is assumed that luminous 
galaxies trace the mass distribution. This means that the parameter b is assumed to 
be unity and, consequently, the relative defects (excesses) of galaxies and dark mat- 
ter are identical inside the void (wall). The case b < 1 has not either theoretical or 
observational support, it corresponds to a void (wall) with an amount of dark matter 
smaller (greater) than that of the case b = 1; by this reason, the anisotropy of the 
case b < 1 is expected to be greater than that estimated in the case 6=1. Similar 
arguments lead to the conclusion that, in the case 6 > 1, the anisotropies are smaller 
than those of the case 6=1. Since the condition 6 = 1 is used along the paper, 
we can state that our computations give upper limits to the anisotropy produced by 
Bootes-like objects, except in the unlikely case 6 < 1. 

In this paper, the wall is described by two parameters, the amplitude, {—) maX i of 
the total density contrast inside the wall and the distance, d w , between the two points 
of the wall in which ^ is the 20% of (—) ma x', the distance d w is called the wall width. 
In the case 6=1, current data suggest that the value of (y) maa; is ~ 4 and the wall 
width d w is ~ Mpc. It can be easily verified (see Section 2) that a wall having 
these features overcompensates the central underdensity described above. The mass 
excess of this wall is about two times the mass defect of the underdense region. This 
is not a model dependent conclusion, but a direct consequence of the observations. In 
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this task, the anisotropy produced by an isolated structure is estimated. The chosen 
overcompensated structure is formed by a void and all the matter surrounding it. 
According to the above observational evidences, overcompensated structures of this 
kind are present in the universe. As required by the cosmological principle, each of 
these structures should be compensated by other structures in large volumes contain- 
ing various voids. In order to imagine this compensation, it is useful to take into 
account that the energy excess surrounding a certain underdensity is shared by the 
neighboring ones and, consequently, this excess also contributes to the compensation 
of other neighboring underdensities; in other words, only a part of this excess must 
compensate the central underdensity. 

The anisotropy produced by a realistic distribution of irregular voids and walls 
cannot be calculated from the anisotropy produced by an isolated void with walls; 
nevertheless, if the chosen isolated structure produces large enough anisotropies and 
its spatial distribution is appropriate, the true distribution of voids could produce 
relevant effects; in this case, the information obtained from the study of isolated 
structures strongly motivates further researches based on suitable approaches. Since 
the effects of isolated voids with walls placed far from the last scattering surface 
are expected to be small, the greatest observed void, namely, the BV has been se- 
lected. Several BV realizations evolving in various backgrounds have been considered 
in order to do an exhaustive study of the anisotropies produced by isolated voids. 
Observational evidences are taken into account in order to select these realizations 
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and backgrounds. 

Henceforth, a is the scale factor, t is the cosmological time, an overdot stands for 
a derivative with respect to t, H is the ratio a/ a, and flo is the density parameter. 
The subscripts D, and B indicate that a quantity has been valued at decoupling, 
at present time and in the background, respectively; for instance, H is the Hubble 
constant. If H is given in units of km s -1 Mpc' 1 , the parameter h = H o /100 is the 
reduced dimensionless Hubble constant. 

The plan of this paper is as follows: Several BV models compatible with current 
observations are defined in Section 2. The initial conditions -at decoupling- lead- 
ing to these models are derived in Section 3. The anisotropies produced by the BV 
models of Section 2 are presented in Section 4; appropriate comparisons with pre- 
vious computations are also given in this section. Finally, the main conclusions are 
summarized and discussed in Section 5. 

2 BV models 

A BV realization is defined by the present density contrast inside the void (^) v , the 
present radius of the underdense region R v , the present amplitude of the density con- 
trast inside the wall {^-) m ax and the present wall width d w . Any present configuration 
of the void in Bootes is a BV realization. Each BV realization is the final state of 
an evolutionary process, which takes place in a certain Friedmann-Robertson- Walker 
Universe. Hereafter, a BV realization and a background define a BV model. Any 
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model describes the time evolution of the void in a certain background. This evolu- 
tion leads to a present state (a realization). If the cosmological constant vanishes, 
h and Q are the free parameters of the background. The numerical codes need a 
fixed value of h; nevertheless, if the distances are given in units of h~ x Mpc, the final 
results do not depend on the chosen value of h; so, only the background parameter 
fl is a physically significant free parameter to be varied. 

Six BV realizations have been selected to be considered in the next sections. 

One of the chosen BV realizations corresponds to (—) v ~ —1, R v ~ 30/i _1 Mpc, 
(^y-)max ~ 2.5 and d w ~ 2h^ x Mpc. This realization is only studied in the case 
r2 = 1- The resulting model is considered with the essential aim of testing our codes 
and comparing our results with previous ones (Panek 1992) 

The remaining five BV realizations are obtained as follows: 

The quantities (^) v and R v are fixed; their values are assumed to be ~ —0.75 and 
30/i" 1 Mpc, respectively. The amplitude (^■) ma x is varied from 2 to 6, and the wall 
width d w is varied from 3k' 1 Mpc to 7hr l Mpc. Taking into account that the values 
of these parameters suggested by the observations (de Lapparent, Geller & Huchra 
1986) are {^) ma x ~ 4 and d w ~ hh^ 1 Mpc, we proceed as follows: in a first step, the 
wall width d w = hh~ l Mpc is fixed and the amplitudes 2, 4 and 6 are considered and, 
in a second step, the amplitude (^) ma x = 4 is fixed and the wall widths are assumed 
to be Mpc, hh~ l Mpc and 7/i _1 Mpc. Note that the realization (^) ma x ~ 4 
and d w ~ hh~ l Mpc is considered in each of the above steps. Hereafter, any of these 
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five realizations is identified by the values of the quantities (— ) max and d w . Each of 
these realizations has been studied in several cases corresponding to tt values ranging 
from 0.2 to 1; nevertheless, for the sake of briefness, only some appropriate models 
corresponding to Q = 0.2 and fl = 1 are presented. 

For each BV realization, the wall compensates the central underdensity at a certain 
distance from the void centre. This distance is called the compensation radius, R c . 
Tables 1 and 2 (seventh column) show these radius for all the models considered in 
this paper. In the cases of the first and fourth rows of Table 1 and the second row 
of Table 2, the compensation occurs outside the wall, while in the remaining cases, 
it occurs inside the wall; for example, in the case (^f-) ma x ~ 4 and d w ~ hh~ l Mpc, 
the compensation radius is R c = 32.7k" 1 Mpc and, consequently, the compensation 
of the central underdensity takes place near the wall centre; this means that about 
one-half of the wall compensates the central void, while the remaining of the wall 
must compensate other underdensities. 

For a given fl value, the quantities defining any BV realization (present contrasts 
and distances) must be numerically obtained -after evolution- from appropriated 
initial conditions; the choice of the initial conditions corresponding to a given model 
(a BV realization evolving in a defined background) is now discussed. 
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3 Initial conditions 

The main goal of this task is the estimation of the secondary gravitational anisotropies 
produced by big voids located far from the last scattering surface. Since the evolu- 
tion of the microwave photons must be studied from the last scattering surface to 
our position in the Universe, the initial conditions for the void evolution are set at 
decoupling time (redshift Zd ec = 1000). 

The initial profiles of the total energy density and the peculiar velocity field fix 
the two arbitrary functions involved in the TBS (see Arnau et al. 1993); hence, these 
initial profiles fix the time evolution of the resulting void from decoupling to present 
time. 

In this paper, the initial density profile is assumed to be 



where R is a radial coordinate and the conditions S\ > 0,^2 < 0, < \s2\ and 
R x i > Ri x are satisfied. 

The initial peculiar velocity is 



1 + 




(1) 



P = P 



1 + {R/R xl f 1 + {R/R x2 f 





where the angular brackets denote a mean value from R = to R. 
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Since the cosmological constant vanishes, the background parameters involved in 
Eqs. (1) and (2) can be written in terms of Q and h. 

Only the choice of the profile (1) is arbitrary. The profile (2) is obtained from Eq. 
(1). It corresponds to vanishing nongrowing modes (Peebles 1980). The form of the 
initial density profile (1) has not any theoretical justification. This form only gives 
a certain parametrization of the initial conditions, this parametrization is expected 
to be suitable in order to describe voids with walls as a result of two facts: (1) 
for small values of R/R x2 , the quantities (R/R x i) 6 and (R/ R x2 ) 6 become very small 
and p becomes quasiconstant; this means that the profile (1) describes a central 
underdensity with a density contrast ~ S\ + e 2 < and, (2) the central underdensity 
is initially surrounded by an overdensity, which is the origin of the present void 
wall. It has been verified that the exponent 6 is suitable in order to get the required 
amplitudes and wall widths at present time, but other exponents could be also tested. 

A realization is defined by the quantities (— )„, R v , (— ) max and d w . The question 
is: which are the values of the parameters £i, e 2 , R x i and R x2 leading to a given 
realization in a fixed background?. 

For each Q value, a numerical code based on the TBS plus Eqs. (1) and (2) 
calculates the quantities (^f-)v, Rv, {^f)max and d w from initial values of ei, e 2 , R x \ 
and R x2 ; this means that the quantities defining a BV realization are not initial 
conditions for the numerical code, but quantities derived from it. 

Given a BV model, the corresponding initial conditions are obtained as follows: 
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arbitrary values of e±, £2, R x i and R x2 are assumed and the resulting values of (—) v , 
R v , (^ £ ) max and dw are compared with those of the chosen model; if they are different, 
the parameters £1, €2, R x i and R X 2 are varied and the results are compared again. 
These calculations and comparisons are carried out by a numerical code based on the 
"gradient method" ( as in Saez, Arnau & Fullana 1993). This code modifies the initial 
values of the parameters £1, e 2 , R x \ and R X 2 in such a way that the new parameters 
lead to a BV model better than the previous one. This code repeats the modification 
of the parameters until the resulting BV model is similar enough to the required one. 
This process requires nonlinear techniques because the final BV model is a nonlinear 
one (see Fig. 1); in other words, nonlinear methods are necessary as a result of our 
BV normalization, which is based on present nonlinear observational data. 

Table 1 gives the initial values of the parameters £1, £ 2 , R x \ and R x2 for each of 
the f2 = 0.2 models studied in this paper. Table 2 gives the same information for 
the Qq = 1 models. The present ^ profiles corresponding to our BV models are 
shown in Fig. 1. For a given realization, the present energy density profiles of the 
models Qq = 0.2 and f2 = 1 are indistinguishable because the same values of (^f)u, 
Rv, (^)max and d w have been chosen in both cases (see Tables 1 and 2). 

It is well known that the TBS only applies before shell crossing. Hellaby and 
Lake (1985) gave the necessary and sufficient conditions for the presence of shell 
crossing; these conditions are satisfied in the cases studied in this paper; hence, the 
shell crossing is unavoidable. The time at which this phenomenon takes place is not 
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given by the Hellaby and Lake conditions; this time must be determined in each 
particular case. In the cases studied in this paper, it has been verified that the shell 
crossing does not take place before present time; hence, the TBS can be used in our 
computations. 

4 Anisotropy 

The initial conditions discussed in Sect. 3 define the space-time structure and, con- 
sequently, these conditions fix the differential equations of the photon trajectories. 
These differential equations must be integrated in order to estimate the anisotropy 
produced by a BV model; this integration is carried out by using the code due to 
Arnau et al. (1993) and Saez, Arnau & Fullana (1993) plus the initial profiles (1) 
and (2). 

The centre of the BV open models (Qq = 0.2) has been located at a significant 
set of distances from the observer, which correspond to redshifts between 0.052 (~ 
150/i -1 Mpc) and 100 (8360/i -1 Mpc); nevertheless, only the results corresponding to 
a few appropriate distances are displayed in the Figures. The centre of the BV flat 
model corresponding to the first row of Table 2 has been placed at the same redshifts, 
while the model used for comparisons with previous computations (second row of 
Table 2) has been only located at lOO/i -1 Mpc in order to facilitate these comparisons; 
hence, all the selected structures are located far from the last scattering surface and, 
consequently, they produce negligible temperature fluctuations and Doppler shifts on 
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this surface; in which, the temperature is assumed to be constant. 

Our code (Arnau et al. 1993, Saez, Arnau & Fullana 1993) numerically computes 
the temperature T of the microwave background as a function of the observation 
angle ip; this is the angle formed by the line of sight and the line joining the observer 
and the inhomogeneity centre. The function T(ip) is then used to calculate the mean 
temperature < T > = (1/2) ^T{ip) simp dip and the total temperature contrast 
S T (ip) = [T(ip) — < T >] / <T>. In the expansion of 8 T in spherical harmonics, 
8t(iP) = D cosip + Q (3cos 2 ip — 1) + higher order multipoles, D and Q 
are the total dipole and quadrupole, respectively. The dipole D is assumed to be a 
Doppler effect appearing as a result of the present peculiar velocity of the observer 
produced by the Bootes-like object; in other words, any gravitational contribution 
to the dipole is neglected; thus the relativistic Doppler quadrupole is D 2 /3. The 
total Doppler effect (dipole and quadrupole) produced by the peculiar motion of 
the observer is subtracted from 5t{iP) to obtain the residual anisotropy 5r{iP) = 
5 T (ip) — D cosip — D 2 /3 (3cos 2 ip — 1); therefore, on account of the large distance 
separating the chosen voids from the last scattering surface, this anisotropy is a 
pure gravitational effect. The rigorous computation of this effect requires nonlinear 
techniques when the amplitude of the density contrast reaches values greater than 
~ 0.1 in some region of the structure. 

The residual anisotropics produced by the BV models of Sect. 2 -for appropriate 
locations- are now presented and discussed. 
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4.1 Open Universe, Qq = 0.2 

The upper panels of Fig. 2 show the residual anisotropy produced by three BV 
models. The background is open (f2 = 0.2), the wall width is d w = 5/i _1 Mpc in 
all the cases, and the values of the amplitude (—) max are 2, 4 and 6. The present 
energy density contrasts of these models are presented in the upper panel of Fig. 
1. The initial values of the free parameters are given in Table 1. In the upper left 
(right) panel of Fig. 2, the void is centred at z — 0.0052 (z = 2.42). The first of 
these redshifts corresponds to the location of the true BV and the second one to 
the position leading to the maximum anisotropy (see below). As it is observed in 
the plot, the greater (^) max , the greater the amplitude of 5 R . In the upper left 
(right) panel, the maximum amplitude of the residual anisotropy is 5 R ~ 8 x 10~ 7 
(5 R ~ 4.5 x 10~ 6 ); this value is obtained in the case (^f-)max = 6. For the model 
d w = 5h~ x Mpc, (^■) m ,ax = 4, the residual anisotropies are 5 R ~ 4 x 10~ 7 (left) and 
5 R ~ 2.3 x 10- 6 (right). 

The bottom panels of Fig. 2 also display the residual anisotropy produced by three 
BV models. The background is the same as in the upper panels, but the realizations 
are different. The wall widths, are 3/i _1 Mpc, 5/i -1 Mpc and 7h~ l Mpc and the 
amplitude is {^ £ -) ma x = 4 in all the cases. The present energy density contrasts of 
these models are given in the intermediate panel of Fig. 1 and the initial conditions 
can be found in Table 1. The redshifts of the void centres are the same as in the top 
panels: z = 0.052 (left) and z = 2.42 (right). In the bottom left (right) panel, the 
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maximum amplitude of the residual anisotropy is 5 R ~ 6.8 x 1CT 7 (5 R ~ 3.8 x 1CT 6 ). 
These values correspond to the maximum wall width d w = 7hr x Mpc. The greater 
d w , the greater the amplitude of 5 R . 

All the cases of the left panels of Fig. 2 correspond to the same backgrounds 
and locations and the central underdense regions have the same structure; hence, the 
differences between the anisotropies of two of these cases are due to the wall. The 
same can be stated for the cases of the right panels. In the cases corresponding to 
the continuous and dashed lines, the anisotropy produced by the wall dominates the 
total effect. 

The upper panel of Fig. 3 shows the residual anisotropy produced by the model 
Qo = 0.2, (^y)max = 4, and d w = hh^ 1 Mpc. Each curve corresponds to a location of 
the void. The redshift defining this location is given inside the panel. This model is 
also considered in Fig. 2. As it is shown in this panel, the amplitude of the residual 
anisotropy is an increasing function of the redshift z -defining the location of the void- 
from z = 0.052 to z ~ 2.42, while it becomes a decreasing function for z > 2.42. The 
maximum amplitude is 2.36 x 10~ 6 , it is is found at z ~ 2.42. As it is pointed out 
below, this behavior is not observed in the case of flat models. 

Fig. 4 shows the density profiles of the three realizations described in the top 
panel of Fig. 1, for z = 2.42 and Q = 0.2. As it can be seen in Fig. 4, these 
structures were evolving in the mildly nonlinear regime when they produced the 
maximum anisotropy. The amplitude of the density contrast inside the underdensity 
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is —0.59 in all the cases. Inside the wall, this amplitude takes on the values 0.48, 0.69 
and 0.82 for the amplitudes 2, 4 and 6, respectively; hence, the standard Eulerian 
linear approach does not apply. 

4.2 Flat Universe, Q = 1 

For Q = 1, the residual anisotropy corresponding to the realization (— ) max = 4, 
d w = hh~ l Mpc is shown in the intermediate panel of Fig. 3. The void centre is 
located at the redshifts displayed inside the panel. The initial conditions are given 
in Table 2. For a flat background, the amplitude of the anisotropy produced by the 
chosen BV realization (the same as in the top panel of Fig. 3) is ~ —3.2 x 10 -7 in 
the case z = 0.052. The modulus of this amplitude decreases as z increases from 
z = 0.052 to z ~ 2.42 and it is an slowly increasing function of z for z > 2.42; hence, 
at redshifts between 1 and 10, the value of this modulus is much smaller than the 
amplitude of the residual anisotropy corresponding to the open case (top panel of 
Fig. 3). At z — 2.42, the ratio between this amplitude and the mentioned modulus 
is ~ 8; therefore, we can state that, at low redshifts, the anisotropies corresponding 
to the open case are much greater than those of the flat case. 

In the intermediate panel of Fig. 3, it can be seen that the wall produces a local 
effect. When the structure is located at z = 0.052, this effect appears between ip = 8° 
and ip = 13°. As the redshift increases, the effect appears at smaller angles. In any 
case, the angular position of the feature coincides with that of the sharped wall. The 
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photons coming along these directions cross a great part of the wall. The accuracy 
of our codes allows us to obtain these small features. 

The sign of 5 R (ip = 0) is positive for f2 = 0.2 and negative for f2 = 1 (see the 
upper and intermediate panels of Fig. 3). In the absence of walls (Arnau et al. 1993) 
as well as in the case of small compensating walls (Panek 1992), the sign of 5 R (ip = 0) 
is negative for both Q values; therefore the sign change only appears in the case of 
overcompensated voids. 

4.3 Comparisons with previous calculations 

Panek (1992) studied four void realizations evolving in a flat background. The model 
(c) of Panek's paper is a BV model having the following features: Q = 1, (—) v ~ —1, 
R v ~ 30k- 1 Mpc, ( 

~~r)max ~ 2.5 and d w ~ 2h ^ Mpc. Our code —based on the 
gradient method- has been used in order to find the initial conditions corresponding 
to this model. The values of e±, e 2 , R x \ and R x2 are given in Table 2 (second row). 
The present density contrast of this model is displayed in the bottom panel of Fig. 
1. It has been verified that the wall compensates the central underdensity at R c = 
51.6/i -1 Mpc. As in Panek's paper, the void centre is placed at lOO/i -1 Mpc from the 
observer in order to compute anisotropies. The residual anisotropy produced by this 
model is plotted in the bottom panel of Fig. 3. This panel and the bottom panel of 
Fig. 1 are to be compared with Figs. 6 and 1 of Panek's paper, respectively. These 
panels have a special format -different from that of the remaining ones- in order to 
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facilitate comparisons with Panek's Figures (Panek 1992). These comparisons clearly 
show that our codes -based on the TBS and Eq. (1) and (2)- have let to a BV model 
very similar to that of Panek (1992); accordingly, the residual anisotropy appears 
to be very similar to that predicted by this author in his case (c). These results 
simultaneously test our codes and those used by Panek (1992). 

Thompson & Vishniac (1987) and Martinez-Gonzalez & Sanz (1990) predicted BV 
anisotropies of a few times ~ 1CT 7 for Bootes-like objects located at z ~ 0.05 in any 
admissible background and similar anisotropies for other redshifts in a flat universe. 
These authors used the Swiss-Cheese model. In the case of small compensating walls 
(bottom panel of Fig. 3) and in the absence of walls (Arnau et al. 1993), our results 
essentially agree with these previous estimates; however, under the following assump- 
tions: (1) overcompensating walls with the features suggested by the observational 
data, (2) an open universe with Q = 0.2, and (3) low redshifts ranging in the interval 
(1,10), previous predictions are magnified by a factor ~ 10. 

5 Conclusions and discussion 

Equation (1) defines a good parametrization of the initial density profiles in the case 
of voids with sharped walls. Our code based on the gradient method -plus Eqs (1) 
and (2)- gives the initial conditions leading to any BV model. 

As a result of the fact that the void walls have been modeled taking into account 
the observational evidences, the compensation of the central underdensity takes place 
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at scales larger than that of a single void. Various voids contribute to this compen- 
sation. 

The anisotropy produced by a Bootes-like void strongly depends on the wall struc- 
ture, the density parameter and the location of the symmetry centre. According to 
previous estimates, which are confirmed in this paper, the anisotropy produced by 
compensated voids evolving in a flat universe has an amplitude of a few times 10~ 7 ; 
however, for f2 = 0.2 and locations between z — 1 and z = 10, the overcompensated 
voids suggested by the observations produce anisotropics of a few times 10~ 6 on scales 
of a few degrees. A question is relevant: What is the effect produced by a distribution 
of voids in an open universe?. In the flat case Q = 1, the anisotropies corresponding 
to the same range of redshifts are much smaller. The value of the density parameter 
is of crucial importance. In the case of overcompensating walls, the sign of the effect 
towards the central region of the void is positive (negative) for open (flat) universes. 
If this effect is detected in future in the case of a single observable structure, results 
could be used in order to constraint the density parameter; nevertheless, it should be 
pointed out that such a detection is not easy, in particular, in the flat case, where the 
resulting anisotropy is very small. 

For the above interval of redshifts (1, 10) and f2 = 0.2, the present distances from 
the void centre to the observer range from ~ 2000/i -1 Mpc to ~ 5900/i -1 Mpc. There 
are many voids located between these distances; nevertheless, only some rare voids 
would be Bootes-like voids (or greater) producing anisotropies of a few times 10~ 6 . 
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Given two observation angles ipi and ip2, the number of big voids n\ and n 2 crossed by 
the photons traveling along the chosen directions can be different. For | rt\ — n 2 |> 1, 
the relative temperature difference corresponding to ipi an d ip2 would be near 10~ 5 . 
This possibility cannot be rejected a priori. It must be either rejected or accepted 
after quantitative calculations. The feasibility of the condition | n\ — n 2 \ > 1 depends 
on the abundance of big voids. Since the anisotropy -on scales of a few degrees- 
observed in experiments as COBE (Smoot et al. 1992) and Tenerife (Watson et al 
1992) is near 10~ 5 , the contribution of big voids at low redshifts could be important. 
In the flat case, this contribution is expected to be too small. In any case, it would 
appear superimposed to the primary anisotropy produced near the last scattering 
surface in the linear regime. 

Similar results were obtained in the case of Great Attractor-like objects (Arnau, 
Fullana & Saez 1994). For Q < 0.4 and 2 < Z < 30, these structures produce 
anisotropies of the order of 10~ 5 on scales of a few degrees. In both cases, either 
the universe is open enough or the anisotropy is negligible. Results about Great 
Attractor-like object enhance the interest of the above question, which should be 
rewritten as follows: What is the anisotropy produced by a distribution of voids and 
great overdensities in an open universe?. 

Although a model of overcompensated isolated voids based on the TBS is currently 
competitive, it has some important limitations related to the spherical symmetry. 
Even if the central underdense region is quasispherical, the true wall is not regular 
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and the motion of the matter contained in this part of the structure is not strictly 
radial. There are clusters and structures in the walls, which produce local peculiar 
motions tangent to the wall. These motions would also produce anisotropy (Tuluie & 
Laguna 1995). The anisotropy produced by the substructures of the walls are expected 
to be important on angular scales smaller than a few degrees and, consequently, the 
estimates of this paper should be admissible. 
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Figure Captions 

Fig. 1. Present density contrast Ap/p(t ) as a function of the present radial distance 
R in units of h~ l Mpc. Upper panel corresponds to three BV realizations with the 
same wall width d w = 5/i~ x Mpc and three different amplitudes displayed inside 
the panel. The BV realizations of the intermediate panel correspond to the fixed 
amplitude (^y-)max = 4 and the wall widths are shown inside the panel in units of 
hr x Mpc. Bottom panel corresponds to (^ £ ) ma x = 2.5, d w = 2h~ 1 Mpc. 

Fig. 2. Left panels show the residual anisotropy 5 R x 10 7 as a function of the 
observation angle ip (in degrees) for several BV realizations placed at z — 0.052. The 
density parameter is flo = 0.2. Upper (bottom) left panel corresponds to the same 
BV realizations as in the upper (intermediate) panel of Fig. 1. Right panels display 
the quantity 5 R x 10 6 for the same models as in the left panels. Void centres are 
located at z — 2.42. 

Fig. 3. Same as Fig. 2. Top panel corresponds to the realization (^■) max = 4, 
d w = 5/i -1 Mpc evolving in an open universe with fl = 0.2. This realization is placed 
at the redshifts displayed inside the panel. In the intermediate panel the realization 
and the redshifts are identical to those of the top panel, but the background is flat. 
The bottom panel corresponds to the energy density profile of the bottom panel of 
Fig. 1. The background is flat. 
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Fig. 4. Density contrast Ap/p as a function of the radial distance R at redshift 
2.42. R is given in units of h^ 1 Mpc. The three BV realizations have the same wall 
width d w = Mpc and three different amplitudes displayed inside the panel. The 
density parameter is fl = 0.2. 
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Table 1. BV models. ^ = 0.2. 



(*£) 

\ p /max 


d w 


£l X 10 3 


e 2 x 10 2 


Rxi x io 2 


i^a x 10 2 


R c 




(h~ l Mpc) 






(h^Mpc) 


(h^Mpc) 


(h' l Mpc) 


2. 


5. 


34.22 


-4.29 


2.60 


2.32 


35.9 


4. 


5. 


5.70 


-1.44 


3.63 


2.09 


32.7 


6. 


5. 


4.91 


-1.36 


4.31 


2.06 


32.4 


4. 


3. 


23.44 


-3.21 


2.70 


2.32 


33.3 


4. 


7. 


4.63 


-1.33 


4.21 


2.05 


32.8 



Table 2. BV models. Sl = 1. 



(*£) 

\ p J max 


dw 


£l x 10 3 


€ 2 X 10 3 


Rxi x 10 2 


i?x2 X 10 2 


R c 




(h~ l Mpc) 






(h^Mpc) 


(h~ l Mpc) 


(h~ l Mpc) 


4. 


5. 


1.43 


-3.63 


3.65 


2.08 


32.7 


2.5 


2. 


0.14 


-12.06 


6.53 


0.99 


51.6 
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